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We study the t-J model with four holes on a 32-site square lattice using exact diagonalization. 
This system corresponds to doping level x = 1/8. At the "realistic" parameter J/t — 0.3, holes in 
the ground state of this system are unbound. They have short range repulsion due to lowering of 
kinetic energy. There is no antiferromagnetic spin order and the electron momentum distribution 
function resembles hole pockets. Furthermore, we show evidence that in case antiferromagnetic 
order exists, holes form d-wave bound pairs and there is mutual repulsion among hole pairs. This 
presumably will occur at low doping level. This scenario is compatible with a checkerboard-type 
charge density state proposed to explain the "1/8 anomaly" in the LSCO family, except that it is 
the ground state only when the system possesses strong antiferromagnetic order. 



It is generally accepted that the physics of high tem- 
perature superconductors is the physics of a doped anti- 
ferromagnetic insulator. Although the detail mechanism 
is still not fully understood, it is believed that frustra- 
tion due to mobile holes may lead to superconductivity 
as well as other competing ordersii^ which may include 
charge ordered and spin ordered phases, or even more 
exotic phases like the staggered- flux phase. One way to 
examine the possibilities of these states is to start from a 
lightly doped antiferromagnet. The t-J model with two 
holes has been extensively studied using various analyt- 
ical and numerical methods 4i^iSiL2i2, These results indi- 
cate that two holes form a weak bound pair with dx^-^y^ 
symmetry. In addition, exact diagonalization study of 
the two-hole model lends support to the staggered-flux 
phase Nevertheless there are still many open questions. 
In a system with two holes the electron momentum dis- 
tribution function (EMDF) (rika) reflects the structure of 
the bound pair only and may not contain information on 
the Fermi surface in the thermodynamic limit. A system 
with four holes is therefore the first step in possibly ex- 
tracting useful information on the system at finite doping 
level. Furthermore, the two-hole system has a few low- 
lying excited states where the holes do not form a dx2_y2 
bound pairi4 It is therefore unclear whether d-wave hole 
pairing still survives when more holes are present. If so, 
how would hole pairs interact? And if not, what scenario 
replaces hole pairing? 

Motivated by the above questions, we numerically di- 
agonalize the t-J model with four holes on a 32-site lat- 
tice with periodic boundary conditions. We caution that 
due to the small size of our system we cannot study long- 
range orders. Therefore superconductivity is not the sub- 
ject of this study. The doping level of this system is 
X — 1/8 which is in the underdoped region. Incidentally 
this is the same doping level as the "1/8 anomaly" in 
the LSCO family where experimental signals for charge 
ordering are more pronouncedt? The Hamiltonian of the 
t-J model is 



n = -tJ2icliij.+h.c.)+j 



TliTlj 



(1) 



where (ij) denotes a nearest neighbor pair. Let us con- 



sider the case J/t = 0.3 first. This parameter is believed 
to be relevant to high Tc cuprates. The ground state has 
energy — 18.682017t, momentum (0,0) and d^^-y^ sym- 
metry. (Note that this refers to the symmetry of the 
total wave function. It is not the pairing symmetry, nor 
does it imply hole pairing.) To see whether there is hole 
binding we calculate the hole-hole correlation functionii 
^(^) — ((1 ~ ^r )(l ~ no)), where is the spin num- 
ber operator in Eq. Fig.Qshows that at J/t — 0.3, 
systems with two and four holes have completely differ- 
ent behaviors. While holes in the two- hole system form a 
bound pair, in the four-hole system they are mutually re- 
pulsive — C (r) increases with r and does not vary much 
beyond r = 2. The root-mean-square separation between 
holes is Tims = (r^) = 2.4732. This is slightly larger 
than 2.3827, which is the root-mean-square separation 
between uncorrelated holes on the same lattice. In order 
to understand this repulsion, we compare C{r) to the pair 
correlation function of four hard core bosons on the same 
lattice. Besides the condition of no double occupancy, 
hard core bosons have effective short-range repulsion due 
to hopping, but are otherwise uncorrelated. As shown in 
Fig. ^ the pair correlation function of hard core bosons 
is strikingly similar to C(r) of our four-hole system. Trms 
of the bosons is 2.4894, which is very close to that of the 
four-hole system. From this similarity we conclude that 
in the t-J model with four holes at J/t = 0.3, holes have 
short-range repulsion which arises from kinetic effect due 
to hopping. 

More information on the interaction among holes 
comes from correlation in momentum space. As usual, 
we define the electron momentum distribution function 
(EMDF) as (nka) = (c^cr'^ko-)- The result is shown in 
Fig. 121 a). Since the electron number is even, (nk|) = 
(nkj^) and we drop the spin variable a from (nko-). In 
interpreting the result it is important to realize that the 
general shape of the EMDF can be misleading. It fea- 
tures a "dome"-shape with a maximum at (0,0) and 
slopes down towards (tTjTt). This is a general feature at 
all doping levels and has been shown^^JS to result from 
minimizing the kinetic energy. It does not represent the 
true Fermi surface. Instead we should focus on those 
k points along the boundary of the antiferromagnetic 
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TABLE L HMDF (rik) in the two- and four-hole systems at 
J /t — 0.3. k points along the AFBZ boundary are marked 
with f . Numbers in brackets are calculated using the additive 
assumption, Eq. ^ and where (rTk)i are from Ref. 
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two-hole system 



0.0031 (0.0058) 
0.0042 (0.0064) 
0.0376 (0.1156) 
0.1234 (0.1416) 
0.1255 (0.1316) 
0.1360 (0.0582) 
0.0497 (0.0137) 
0.0104 (0.0021) 
0.0657 (0.0169) 



four-hole system 



0.0041 (0.0115) 
0.0103 (0.0128) 
0.1813 (0.2311) 
0.2565 (0.2832) 
0.2529 (0.2633) 
0.2454 (0.1164) 
0.0594 (0.0274) 
0.0185 (0.0042) 
0.0970 (0.0339) 



FIG. 1; (Color online) Hole-hole correlation in the ground 
states with two and four holes. The thick shaded line is the 
pair correlation function of four hard core bosons on the same 
lattice. 
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FIG. 2: EMDF (nk> in the ground state of the four-hole 
system at (a) J/t = 0.3 and (b) J/t = 0.8. Solid line is the 
AFBZ boundary. Due to symmetry, only one quadrant of the 
Brillouin zone is shown. 



Brillouin zone (AFBZ), i.e., from (tt, 0) to (0,7r). They 
are not affected by the kinematic effect and therefore 
reflect the physics of the system.'^ Fig. shows that 
along the AFBZ boundary the EMDF has a minimum at 
(7r/2,7r/2). Again this is in sharp contrast with the two- 
hole system whose EMDF has a maximum at (7r/2,7r/2) 
along the AFBZ boundary.^ Unfortunately we do not 
have useful information on how the EMDF changes from 
(0,0) to (7r,7r) because kinematic effect dominates along 
this direction. Therefore we cannot conclude whether 
the minimum at (7r/2,7r/2) implies a hole pocket or a 
half pocket (segmented Fermi surface). But it is not in 
contradiction with either one. 

The fact that the EMDF of the four-hole system has 
a minimum at the single-hole ground state momentum 
(7r/2,7r/2) immediately suggests the relevance of the 
single-hole ground state to the multiple-hole oneii^si^ii^ 
If holes doped into the parent system behave like weakly 
interacting fermions, then it is reasonable to expect that 



many-hole systems can be approximated by filling up the 
single-hole band. This should lead to "hole-pockets" at 
single-hole ground state momentafiS We would like to see 
how well the EMDF of the four-hole system agrees with 
this simple picture. For this purpose we consider the hole 
momentum distribution function (HMDF), 



(nicer) 



{N„ + Nh)/N- 



(2) 



where iVo- and Nh are numbers of a spins and holes re- 
spectively. If multiple-hole states can be built up from 
single-hole states, we expect their HMDFs to be additive, 



(nk)2 ^ (nkT)i + (nkx)i, 
(nk)4 ^ 2((7ikT)i + (?7ki)i). 



(3) 
(4) 



where (rlkcr) jv,, is the HMDF of the Nh-^o\e system. Ta- 
ble shows the HMDF in the two- and four-hole systems 
together with the approximations from Eq. ||2J) and Q 
respectivelyjii Note that Eq. (O and Q lead to very 
prominent maxima at (7r/2, 7r/2) along the AFBZ bound- 
ary. This of course does not agree with our exact re- 
sult in the two-hole system. Even in the four-hole sys- 
tem, the maximum is too sharp compared to the exact 
result. Therefore a simple additive model as given in 
Eq. © and 10} does not work in the t-J model. We 
remark that such simple additive picture holds in the ex- 
tended t-J model with parameters chosen to represent 
electron-doped cuprates>i^ We conclude that although 
the EMDF of the four-hole system exhibit minima that 
suggests the relevance of the single-hole ground state to 
the multiple-hole one, there must still be appreciable cor- 
relation among the holes. 

Next we look at a subtle form of correlation that 
is believed to reflect hole binding. A staggered pat- 
tern in the current correlation function {jkijmn), where 
jki — ii(c^Ci — cjck) is the hole current on the nearest 
neighbor bond kl, was first suggested as a manifestation 
of the staggered flux phasepiS, Although there exists an- 
other explanation in terms of the spin polaron theoryj^S 
both theories point to the fact that such staggered cor- 
relation reflects hole binding in the d-wave channel. The 



3 




FIG. 3: Current correlation in the ground state of the four- 
hole system at J jt — 0.3. The reference bond mn is indi- 
cated by a shaded line. Except for the reference bond, line 
widths are proportional to (jkijmn) /x, where x = 1/8 is the 
doping level. For reference purpose, the numerical values of 
(jkijmn) /x are shown next to two horizontal bonds. 



current correlation at J/t = 0.3 is shown in Fig. |3| We 
note that the direction of the correlation function on 
many bonds deviate from a staggered pattern. It cer- 
tainly does not have the clear staggered pattern that ex- 
ists in the two-hole systemiiS Another way to display 
the overall pattern is through the vorticity correlation, 
Cvvir) = {V{r)V{0))/x. The vorticity V{r) of a square 
plaquette centered at r is defined by summing up the 
current correlations around it in the counterclockwise 
direction.^- We plot Cvv{r) in Fig. ^ It clearly shows 
that the vorticity correlation is weak and decays rapidly, 
leading to deviation from a staggered pattern at large dis- 
tance. Hence from the current correlation we conclude 
that there is no sign of d-wave pairing in the four-hole 
system at J/t = 0.3. 

Up to now we have shown that when the number of 
holes in a 32-site lattice is increased from two to four, d- 
wave pairing no longer exists. We note that at the same 
time antiferromagnetic (AF) correlation in the spin back- 
ground is also lost. Evidence comes from Fig. O which 
shows the spin correlation function (So • Sr) and static 
structure factor S'(k). At J/t — 0.3, when we increase the 
number of holes from two to four, the spin correlation de- 
creases very rapidly and has large fluctuation. The sign of 
the correlation at larger distances shows deviation from 
Neel order. In the momentum space this leads to a strong 
suppression of the AF peak in S'(k) at (tTjTt), leaving no 
sign of AF correlation at all. Note that this is consistent 
with the phase diagram of La2_xSra;Cu04, where a small 
doping level of a few percent is enough to destroy the AF 
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FIG. 4: The vorticity correlation Cvvir) in the ground state 
of the four-hole system at J/t = 0.3. The reference plaquette 
at r = is indicated by a cross inside it. Because of symmetry, 
only a quarter of the cluster is shown. 

phase. This partially answer the question we posed at 
the end of the first paragraph in this paper — when the 
doping level is large enough to destroy AF correlation, 
holes do not form bound pairs. Instead they become mu- 
tually repulsive. Then how about at lower doping level 
where AF correlation still exits? Do our two-hole results 
really represent the behavior of a lightly doped system, 
or are they just due to the artifact that there is only 
one hole pair? To answer this question directly one has 
to solve the t-J model with a few hole pairs on a larger 
lattice. But this is currently not feasible. In order to 
mimic a lightly doped system with a few hole pairs, we 
put two hole pairs on a 32-site lattice with an AF spin 
background. For this purpose we increase J/t to 0.8. The 
ground state at this parameter has energy — 31.917820t, 
momentum (0,0) and s symmetry. Fig. El^a) shows that 
its spin correlation decays almost as a power law, show- 
ing that it has at least local antiferromagnetic order. Its 
static structure factor clearly has a peak at (tt, tt). We 
therefore believe that the four-hole system at J/t — 0.8 
possesses AF order. Next we are going to investigate the 
possibility of d-wave pairing in this system. 

A first sign of d-wave pairing at J/t = 0.8 comes from 
the current correlation shown in Fig. El Except for a 
few bonds at or near the edge of the lattice, the direc- 
tions of the current correlation at all other bonds show 
the correct staggered pattern. Compared to the result 
at J/t — 0.3 (Fig. EJ, the current correlation is much 
stronger in the present case, especially at small distances. 
When we sum up the current correlation around plaque- 
ttes to obtain the vorticity correlation, a clear staggered 
pattern emerges. Fig. [3 clearly shows that the vortic- 
ity correlation has a staggered pattern with a charac- 
teristic TT phase shiftfi^ i.e., the sign of {V{r)V{0)) is 
(^_lYx+ry+i _ Note that the vorticity correlation decays 
quit fast, which is a characteristic of strong hole binding 
at J/t = 0.8ii£ This shows that the four-hole system at 



4 




FIG. 5; (Color online) (a) Spin correlation function (So ■ 
Sr), and (b) static structure factor 5'(k). Empty and filled 
symbols in (a) represent positive and negative correlations 
respectively. 



J /t — 0.8 possesses at least short-range staggered current 
correlation. A second sign of d-wave pairing at J jt — 0.8 
comes from the hole correlation C(r) in Fig. It con- 
sists of two regions with different behavior. At r < 2 it 
resembles that of a two-hole system which shows a de- 
caying trend and with a maximum at r = V2. The latter 
feature has been shown to be due to the d-wave nature of 
the two-hole bound stateiSLSS At r > 2, C(r) remains al- 
most constant, showing that holes at larger distances are 
uncorrelated or even repulsive. This suggests that holes 
form d-wave bound pairs, and the two hole pairs repel 
each other. Further evidence for this scenario comes from 
the EMDF. First of all we note the qualitative similar- 
ity between the EMDFs of the four-hole [Fig. Efb)] and 
two-hole systems (Fig. 8 of Ref. 0) — along the AFBZ 
boundary they are largest at (7r/2,7r/2) and smallest at 
(37r/4, 7r/4). If the two hole pairs like to stay away from 
each other, then their mutual influence will be a mini- 
mum. Consequently we expect that the EMDF of a sys- 
tem with two hole pairs can be built up from that with 
one pair in an additive manner. This is demonstrated 
in Table HH which shows that the HMDF of the four-hole 
system along the AFBZ boundary is roughly twice that 
of the two-hole system. 

To summarize, our results indicate a close relation be- 
tween d-wave pairing and AF correlation. They favor 
the scenario that holes in the lightly doped t-J model 
form weak d-wave bound pairs when there is AF corre- 
lation in the spin background, and there is an effective 
repulsion among the d-wave pairs. In this case the sys- 
tem has all the characteristics previously reported for a 
system with two hole3^i^ii& — the hole correlation func- 
tion has a maximum at r = 72, and the EMDF IS a 
maximum at (7r/2,7r/2) along the AFBZ boundary. The 
current correlation function exhibits a staggered pattern. 
All these characteristics reflect the d-wave nature of the 




FIG. 6: Same as Fig. IHlexcept at J/t = 0.8. 
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FIG. 7: Same as Fig.glexcept at J/t = 0.8. 



TABLE II: HMDF in systems with two holes ((nk)2) and 
four holes ((nk)4) at J jt — 0.8. k points along the AFBZ 
boundary are marked with f . 
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bound state. At higher doping level where AF correla- 
tion disappears, holes no longer form bound pairs and 
they become mutually repulsive. This repulsive behav- 
ior is very similar to that of hard core bosons where the 
repulsion is short range and is due to lowering of kinetic 
energy. The EMDF exhibits dips at (7r/2,7r/2) along 
the AFBZ boundary. This agrees with mean-field theory 
predictions, although we cannot determine whether 
the dips are hole pockets or segmented Fermi surfaces. 
We emphasize once more that since we cannot study long- 
range correlations in our small system, we do not know 
whether those d-wave hole pairs can lead to superconduc- 
tivity. Therefore although our result indicates the coex- 
istence of AF and d-wave pairing, it does not necessarily 
imply the coexistence of AF and superconductivity. 

One popular candidate of a charge-ordered phase in a 
doped antiferroniagnet is the stripe phase where charge 
carriers form one-dimensional strips. It is used to ex- 
plain the "1/8 anomaly" which, incidentally, is at the 
same doping level as our four-hole system. But since the 
model we use is isotropic in the x-y plane, we are not 
able to detect any anisotropic correlation. Therefore our 
results are inconclusive regarding stripe phase. Another 
possible charge-ordered phase that can give rise to en- 
hanced charge order at x = 1/8 is the checkerboard- type 
ordering of hole pairsiSSiS& This ordering emphasizes the 
two-dimensional nature of the spin system in contrary 
to the one-dimensional nature of the stripe phase. In 
the checkerboard-type picture, holes form d-wave bound 
pairs which in turns form charge-ordered states at a series 
of magic doping levels. At x = 1/8 these bound pairs are 
arranged in a square lattice of sides 4. Incidentally this 
configuration is compatible with our results when there 



is AF correlation — the long distance correlation of hole 
pairs in the checkerboard-type order will make hole pairs 
appear to be repulsive on a small lattice. As a result 
the hole correlation will have characteristics of d-wave 
pairs at small distance and becomes repulsive at larger 
distance, very much like the one at J /t = 0.8 in Fig. ^ 
But according to our results the system must possess at 
least local AF order in order for this to be possible. At 
a; = 1/8 and realistic value of J/t, AF correlation no 
longer exists and such a charge-ordered phase should not 
be the ground state. Nevertheless, this charge-ordered 
phase may exist as an excited state. Therefore our re- 
sults can be interpreted as providing indirect evidence 
for the checkerboard-type charge ordering. But we cau- 
tion that this does not mean that our results favor the 
checkerboard-type ordering to stripe phase. Due to its 
isotropic nature, the model we use naturally favors two- 
dimensional ordering to one-dimensional one unless other 
effects (such as lattice distortion) play a role in selecting 
the ground state. Finally we remark that our results do 
not seem to support phase separation at this doping level 
because we do not see evidence that four holes tend to 
cluster together even when the spin exchange integral J 
is as large as 0.8i. 
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Note added - after this work is finished we learn that 
A. Laeuchli has independently solved the t-J model with 
four holes on a 32-site lattice. 
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